New 1/N expansions in random tensor models 

Valentin Bonzom 1 ^! 

1 Perimeter Institute for Theoretical Physics, 31 Caroline St. N, ON N2L 2Y5, Waterloo, Canada 

(Dated: November 8, 2012) 

Although random tensor models were introduced twenty years ago, it is only in 2011 that Gurau 
proved the existence of a 1/N expansion. Here we show that there actually is more than a single 
1/iV expansion, depending on the dimension. In the large N limit, these new expansions retain 
more than the melonic graphs. Still, in most cases, the large N limit is found to be Gaussian, and 
therefore extends the scope of the universality theorem for large random tensors. Nevertheless, a 
scaling which leads to non-Gaussian large N limits, in even dimensions, is identified for the first 
time. 
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o 

> \ INTRODUCTION 

o 

Random tensor models were introduced twenty years ago as a generalization of random matrix models [IH3 ■ While 
£ — . ' a matrix is a two-dimensional array, a tensor is a d-dimensional array, d > 2. Partition functions of tensor models have 
an expansion onto Feynman graphs which arc dual to triangulated <i-dimcnsional pseudo-manifolds, and arc therefore 
of interest for the path integral quantization of general relativity. They indeed realize a sum over discrete geometries 
and the hope is to define quantum gravity by taking the continuum limit. 

However, the 1/N expansion of matrix models, which organizes the sum over two-dimensional topologies at large 
matrix size (i.e. small Newton constant) has been missing in tensor models, until Gurau eventually found one 
i in 2011 |5|— lV| . This 1/N expansion applies in particular to the set of complex tensors with measures of the form 

exp — S where S is a combination of tensor invariants with some U(N) symmetry properties It has triggered a 
full analysis of tensor models, which reproduces the key features of matrix models. The dominant contributions at 
large N have been identified (some triangulations of the d-sphere) and proved to be summable Their singular 
\ behavior provides a continuum limit and coupling to (non-unitary) matter generates different universality classes in 
■ the continuum whose critical exponents can be found exactly [1(1 EH ■ 

' While these developments arc analogous to those obtained in matrix models, and therefore put tensor models 
to a somewhat similar level for the first time, the properties of the large N limit are quite different. First, the 
multi-critical behaviors in the continuum limit correspond to the universality classes of random branched polymers 
(entropy exponents 7 = 1 — 1 /m in contrast to 7 = — 1/m in matrix models). This was in fact expected from numerical 
simulations, performed in the time when no 1 /N expansion was known in tensor models. In those numerical simulations 
on triangulated spheres in dimensions d > 3 (a program called Euclidean Dynamical Triangulations) , it was observed 
£> ■ that at large N (small Newton constant) the partition function is dominated by stacked spheres, whose statistical 
J^j ! properties are those of random branched polymers flil - fbsj ]. This has been confirmed analytically using the large N 
. limit of tensor models [9l-ITlj . 

Second, random tensor models become Gaussian at large N when d > 3. This is a strong universality theorem, 
initially proved in [I2J (then in [l]| by a different method), which does not hold for random matrices. While it does 
not prevent critical behaviors in the continuum limit (the full large N covariance crucially depends on the details of 
the model), this theorem directly leads to the critical exponents of random branched polymers. 

The origin of the universality theorem is to be found in the structure of the dominant Feynman graphs at large N. 
It only probes the melonic sector, where the Feynman graphs dominating the free energy maximize the number of 
faces (for the triangulated sphere it means maximizing the number of (d— 2)-dimcnsional simplices at fixed number of 
top simplices), and the melonic Feynman graphs actually create more faces than the planar ribbon graphs of matrix 
models (at fixed number of vertices and lines). This melonic sector has the same universality class as the stacked 
spheres [9], the one of random branched polymers. Moreover their combinatorics is so specific that only Gaussian 
behaviors can survive at large N. 

Therefore, both the statistical properties in the continuum and the universality theorem are tied to the fact that 
the large N limit only probes a very specific subset of triangulations. However, there is a priori more freedom in 
tensor models than in matrix models. The reason is that in matrix models, there exists only one invariant at order 2p 
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in the matrix elements, ti(MM^) p , whereas many more exist when d > 3 Q. (Physically, those invariants represent 
boundary triangulations. In 2d, they can only be loops, but as the dimension goes to three and higher, more boundary 
triangulations exist.) The question is thus whether other 1/N expansions exist in tensor models which would probe 
more invariants than the melonic ones. 

• The first result of the present article is that random tensor models admit new 1/N expansions. They are based 
on the same ingredients as Gurau's, but in the large N limit, they retain more than the melonic graphs. 

• However, in most cases, the large N limit is still Gaussian. Although it might look surprising at first sight, the 
reason is that dominant contributions must have melonic subgraphs and the universality result applies to them. 
This is the second result, which extends the scope of the universality theorem. 

• The third result is that in even dimension, there exists a 1/N expansion which is not Gaussian at large N, 
because it is based on planar instead of melonic subgraphs. We prove it exists but do not solve this case in the 
present paper and leave it for future investigations. 

The article is organized as follows. In the Section U we review Gurau's 1/N expansion, and the solution at large 
N. The same notions are then used in the Section |n] to make sense of new 1/N expansions. Their large ./V limits are 
discussed in the Section IlIIl including the extension of the universality theorem as well as the new large N limit in 
even dimensions which is not Gaussian at large N. 



I. GURAU'S 1/N EXPANSION 

The presentation follows Q, to which the reader is referred for more details and references. 



A. Invariants, colored graphs and their degree 

A tensor T with components T ai ... ad , a% = 1, . . . , N, is a multi-dimensional array of N d complex numbers. To study 
large random tensors, we need a family of invariant objects constructed from its components. Let us consider the 
action of the external tensor product of U(N), 

T ai ---a d = ^aifci " ' U ad b d — &d! (1) 

bi,...,b d 

for a set ofd independent unitary matrices {U^}, while T transforms similarly with {j/W}. The use of different copies 
of U (N) acting independently on the different indices implies that invariant quantities are obtained by contracting 
an index of a T in position i with an index of a T in the same position i. Taking p tensors T and p tensors T, we 
saturate all indices two by two, respecting their positions. Products and sums of such invariants generate the full 
algebra of invariants. 

It is standard to represent them graphically in terms of line-colored graphs called bubbles. To each T, we associate a 
white vertex and to each T a black vertex. The indices are represented as d half-lines attached to the vertex. Because 
the position of the indices must be respected in the contractions, half-lines carry a label which is just the position 
of the corresponding index. This index is called the color of the line and can take the values i = 1, . . . , d. When an 
index is contracted between a T and a T, the two corresponding half-lines (which have the same color) are connected 
to form a line between a white and a black vertex. 

When all indices are saturated, a closed, bipartite graph with vertices of degree d and colored lines such that each 
vertex has all d colors on its lines, is obtained. Examples are given in the Figure [T] The connected ones generate 
the algebra of invariants. In the remaining, bubble refers to such a bipartite, connected, d-colored graph, and the 
corresponding tensor invariant is denoted B(T,T). 

At d = 2 there is a single bubble with 2p vertices, the loop with alternating colors 1,2, and it corresponds to the 
invariant tr(TT+)*\ 

Let / be a finite set and {-Bj}ie/ a finite set of bubbles which we use to define an invariant action 



S(T,T) = T-T + ^2t i B i (T,T), 
iei 



(2) 



Figure 1. Tensor invariants can be represented as bubbles, i.e. closed bipartite d-colored graphs. Here d = 3. 




Figure 2. In a Feynman graph, the vertices of the bubbles receive an additional line which connects a white vertex to a black 
vertex. They are represented as dashed lines with the color 0. Therefore Feynman graphs are (d + l)-colored graphs. Here 
d = 3. 



where {U} is a set of coupling constants, and T ■ T is the natural contraction of T with T, corresponding to the bubble 
with two vertices and d lines connecting them. The partition function and the free energy are 

Z{{ti}) = cxp -F({ti}) = J [dT df] exp -N^ 1 S . (3) 

The free energy F admits an expansion onto closed connected Feynman graphs built from the bubbles Bi . Each vertex 
of a bubble Bi is a T or a T, and therefore, due to Wick's theorem, they receive an additional line, the propagator. 
Each of these new lines connects a T to a T. By giving them the fictitious color 0, it turns out that the connected 
Feynman graphs which enter the expansion of the free energy are just (d + l)-colored graphs, as can be seen in the 
Figured 

The amplitude that each graph G receives has a scaling with TV which is the key interest of this paper. The 
propagator is 

C ai ---a d .b 1 ---b d ~ jyd-1 ^<»1&1 ^a d b d - (4) 

Therefore, each line of color brings a factor N~ d+1 and the number of such lines is the half-number of vertices of G, 
denoted p. Each bubble brings a factor N . Moreover, the Kronecker deltas of the propagator and the contraction 
of indices along the colored lines of the bubbles identify the indices all along the loops with alternating colors and 
/x, for [A = 1, . . . , d. We call these loops faces with colors 0/x and denote /o M their number. There is a free sum over 
= 1, . . . , N for each of them, so the amplitude of a connected graph G with 2p vertices is 

A{G) = N^Ufo.-(d-i)p+(<i-i)b Y[t1. (5) 

iei 

Here bi denotes the number of bubbles Bi in G, and b = ^2iei ^» ^ ne total number of bubbles. The key ingredient to 
control this scaling is a formula which counts the number of faces in terms of the number of vertices in any connected 
colored graph. For ribbon graphs, we know that there is a topological balance between faces and vertices (and lines) 
encoded in the genus of the ribbon graph Here the structure of colored graphs is probed by identifying some 
ribbon sub-graphs called jackets. A jacket goes along all vertices and lines of G but only a subset of faces (fa-El • 

Definition 1. Let A > 3 be an integer, G be a connected, A-colored graph with 2p vertices (hence Ap lines) and a 
be a cycle on {1, . . . , A}. We call the loops with alternating colors \x, v , for fj, ^ v both in {1, . . . , A}, faces of colors 
/i, v. The jacket J associated to a is the connected ribbon graph which contains all the faces of colors (cr 9 (l), o~ q+1 (Vj) 
for q = 0, . . . , A — 1 in G. Therefore the number of faces in J is given by fj = 2 — 2gj + Ap — 2p, where gj is the 
genus of J. We define the degree U)(G) of G as the sum of the genera of the jackets. 
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Notice that the degree is a positive number. Since there are |(A — 1)! distinct jackets and every face appears in 
exactly 2(A — 2)! jackets, a sum over all jackets allows an over-counting of faces, leading 0-0] to the following 

Theorem 1. For G a connected A-colored graph with 2p vertices and fa faces 

fe= (A- 1 )(A-2) p+(A _ 1) __^_^ (c) {tj 

Therefore, for our connected (d+ l)-colored graph G, 

E/o„ = / G -£ Mb* = ^J^P + d -uTYV^G) J2 (7) 

When d = 2, one simply has /b ; = 1 and the 1/N expansion of matrix models is reproduced straightforwardly. When 
d > 3, the bubbles are connected d-colored graphs, so that the formula ([6]) applies to the number of faces , leading 
to 

A, d(d-l) , 2 ,^ , r(d-l)(d-2) , 2 

E/0M = -4- J p + d-( S ri5 i «(G)-E^ i f ^ + £i - 1 -(^^ 

= (d - l)(p - 6) + d - ^r^(G) + ^jj E 6i 

It is shown in [ijj that ^J^i ^G) — rgjfjl Siei ^* ^(Bi) i s a positive integer. As a result, the free energy admits the 
1/N expansion 

F({U}) = N d V iV-T^^+T^E,,^,^) TT^ (g) 

— ' s(G J 

connected graphs G ££i 
with d + 1 colors 

where s(G) is a symmetry factor. 

Similarly, the expectation values of bubbles, defined by 

(B(T,T)) =e F J [dTdT] B(T,T) exp —N d ~ 1 S(T, T), (10) 

admit an expansion onto connected, (d + l)-colorcd graphs which contain B as a marked subgraph. The scaling of 
each graph is obtained by assigning the scaling derived above for a (d+ l)-colorcd graph and remove the contribution 
which comes from the insertion of B as a bubble, i.e. N^ 1 . Thus, 

connected graphs G i&I 
with d + 1 colors 
containing B 



B. The large TV limit: melonicity and universality 

The scaling entering ([9]) can be re-expressed as 

(d^Tj! u & - Jdhy. g ^ ^ = l^ G) + dtdhy. ( W(G) - rf g ^ ' (12) 

where w(G) — d^,- g/ biUi(Bi) can be proved to be positive [ItJ (and proposition 2 in Appendix A of Q). Therefore 
the scaling is a sum of positive terms, and at large N only the graphs for which they vanish can survive. As a result , 
the free energy in the large N limit is dominated by graphs G with vanishing degree, built from bubbles Bi whose 
degrees are also zcrcP|. This is in contrast with matrix models in which all invariants tr(TT^) p can contribute to the 
large N limit. Here only the bubbles with ui(Bi) = are relevant in the action. 



1 Note that u>(G) = implies that the degrees of its bubbles vanish. 



(a)An elementary melon (b)A closed melonic bubble with 4 colors, 

with external color 1. It 
has 2 vertices and creates 
A(A - l)/2 new faces. 

Figure 3. Melonic colored graphs are built by recursive insertions of the elementary melon. 




Figure 4. Part of a melonic bubble (solid lines), at d = 4, where (v, v) is a canonical pair. It means that the insertions Ii,l2, 13 
are melonic 2-point insertions. All melonic graphs with the additional color (the propagator, drawn in dashed lines) are 
obtained by connecting the vertices of the canonical pair with any melonic 2-point insertion Iq. 



Bubbles and graphs of degree zero are called melonic and built as follows. From the formula © it is found that 
lu = corresponds to maximizing the number of faces at fixed number of vertices. In addition, to keep the degree 
constant when adding two vertices to a graph, A(A — l)/2 faces have to be created. This can be realized locally by 



the 2-point insertion pictured in the Figure 3(a) where all lines of the vertices are paired but the two external ones 
(which have the same color). It is called the elementary melon. Then it can be proved that all melonic graphs are 
obtained by recursive insertions of this patch on any line, starting from the closed graph with two vertices An 
example is given in the Figure |3(b)| 

In the large N limit of the model, only melonic bubbles survive and the (d + l)-colorcd graphs G built from them 
must be melonic too. Notice that from the way melonic colored graphs are built, the vertices are inserted as canonical 
pairs. To identify them within a given melonic graph, one proceeds by recursive removals of the elementary melons. 
The graph contains some elementary melons with obvious canonical pairs. To identify the others, one removes the 
elementary melons and replaces them with lines. This gives another melonic graph on which the canonical pairs of 
the elementary melons arc trivially identified, and so on, until one is left with the graph with exactly two vertices. 

Given a melonic bubble B with 2p vertices, there is only one way to build melonic (d + l)-colorcd graphs using 
propagators and the bubbles {Bi} of the action, which is by connecting the vertices of the canonical pairs with lines of 
color and inserting any melonic 2-point function on them The Figure 2] illustrates these melonic contractions. 

The full expectation value of B(T,T) is thus obtained by inserting the full 2-point function (T ■ T) between the p 
canonical pairs. This explains the large N universal behavior, 

(B(T,f)) = N [U({ti}] p , with U({ti}) = (T-T)/N, (13) 

for melonic bubbles with 2p vertices, and lirn/v'_ >00 (.B)/./V' = otherwise. Therefore these models are Gaussian at 
large N. To complete the solution, the completely dressed covariance U has to be found as a function of the coupling 
constants. This is done using a Schwinger-Dyson equation, which turns into a polynomial equation Q, 

1-17-5^^^ = 0, (14) 



where pi is the half- number of vertices of B{, i G /. 
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II. THE NEW 1/N EXPANSIONS 

a. Color slices. In a (d + l)-colored graph, the subgraphs formed by all the lines with colors in a subset of 
{0, ...,d} are also colored graphs and therefore have a degree. The idea of our new 1/N expansions is to make the 
Feynman graphs scale with the sum of the degrees of such colored subgraphs, instead of the degree of the full graph. 
To realize this, we need to chop the set of colors (index positions in the tensor) {1, . . . , d} into L subsets denoted Di, 
I = 1, . . . , L. Since it does not matter which color is in which subset (only the number of colors will matter), we keep 
the natural order on {1, . . . , d} and cut that set into slices. Choose a set of integers {c?;}; = o....,l such that d;_i < di 
and do = 0, dx = d. The color slice D[ is the set containing all colors between d;_i + 1 and di, 

Di = {d;_i + l,...,d,}, 1 = 1,. ..,L. (15) 

It contains A; = di — df_i colors, for I = 1, . . . , L. The set of colors is partitioned into these L slices, 

L 

{l,...,d}=\jD l . (16) 
i=i 

Each color slice contains at least two colors. Note that at fixed d, L can take values between 1 (for which D\ contains 
all colors), and d/2 for d even (all Di have exactly two colors), and \_d/2\ + 1 for d odd (all Di have exactly two colors 
except one with three colors). 

At d = 2 and d = 3, this is trivial (since color singlets are not allowed): L = 1 and there is only one subset which 
contains all the colors. It becomes interesting for d > 4. When d = 4, we have two possibilities, D\ = {1, 2, 3, 4} with 
L = 1, but also D\ = {1,2},L>2 = {3,4} with L = 2, up to color relabeling. Similarly for d = 5, Di = {1,2,3,4,5} 
with L = 1, and D\ = {1, 2}, D2 = {3, 4, 5} with L = 2 up to color relabeling. 

b. Classification of colored graphs. This slicing induces a natural classification of bubbles (connected d-colored 
graphs) and Feynman graphs (connected (d + l)-colorcd graphs) according to their number of connected components 
in the slice For a bubble B, the subgraph made by all the vertices and all the lines with colors in D; is denoted 
B"' and its number of connected components n(B^). Each connected component is a connected A;-colored graph. 

For a connected (d + l)-colored graph G with colors 0, . . . , d, G"> is defined as the (A; + l)-colored subgraph made 
by all the vertices and all the lines with colors in D\ plus the color 0. It has n(G">) connected components. 



A. Introduction of scalings on the slices 

Now we choose the scaling of our tensor models as Gurau's scaling in each color slice. We know that for connected 
bubble invariants with A colors, the Feynman graphs get A + 1 colors and the correct scaling in front of the action is 
jyA-i_ k ere j s |- apply this principle to all color slices Di, each contributing N Al ~ 1 , and superimpose their 

scalings, yielding 

L 

j~J jyA,-l _ x jyrfa-di-l x . . . x jyd-d L -i-l _ _/V d ~ L . (17) 

1=1 

Let / be a finite set. We consider a finite set of connected, d-colored bubbles {i?i} !e / such that they have only one 

connected component in the slice i.e. n (Bf ] ) = 1. Remarkably we will be able to relax this requirement in the 
Section [II Bl but keeping it for the time being makes the presentation simpler. The action is 

S = T-f + J2 t i B i( T ,T), (18) 
iei 

which has no explicit TV-dependence. The partition function and the free energy are 

Z({U}) = exp-F({U}) = J[dTdT] cxp-N d - L S . (19) 

For L = 1, this reduces to the scaling of the Section fl] 

Let us show that the above scaling leads to a 1/N expansion. The Feynman graphs entering the expansion of the 
free energy are connected (d + l)-colored graphs built from the bubbles Bi. The amplitude A(G) of a graph G with 
2p vertices is found by changing d — 1 to d — L in the Equation ([5]) , 



(20) 
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The number of faces needs to be re-expressed in terms of the number of vertices. This is done independently in 
each color slice. Since the bubbles Bi have only one connected component B^' in each slice and G is connected, the 
subgraphs are also connected. Therefore the counting of faces in the slice / is performed by applying the Equation 



where 



£ fo, = (A, - l)(p — &) + A; - ^-^( G CD) + j^—- £ h U (B®), (21) 



^-fATwEMBf)^. (22) 



A(G) = 7V d -^=^T^"( G<!> ) + T^^^^( B . (i) )J TT t J* ) (24) 



(Ai-1)! v 7 (Ai-2)!^ 
Summing over the slices I = 1, . . . , L, we get 

^/o, = (d-i)(p-6)+d-Ef7A^ w(G(0) + 7A^^ 6 ^ ( ^ 0) ■ 
ju=l i=i L \ 1 >' \ 1 r i £ j 

Therefore the amplitude of a graph is bounded by N d , 

and the free energy has the following 1/N expansion 

F({U}) = N d J2 N-^TK^^^+^^^b^iB^)] (zi^-Q^. ( 25) 

connected graphs G ^ 
with d + 1 colors 

B. Extension to arbitrary invariants 

In the above section, only bubbles with exactly one connected component per slice I were considered. We now relax 
this requirement and show that our new scalings allow for arbitrary bubbles. The danger is that some bubbles, like 
the melonic ones, create too many faces at fixed number of vertices with respect to our scaling so including them 
would make the free energy unbounded. However, a well-chosen, bubble-dependent, large -/V suppression can be added 
in the action to make the graph amplitudes scale uniformly. 

The bubbles Bi, i G /, can now have an arbitrary number n(B^) of connected components in the slice I, each 
being a A;-colored graph. The degree u (Bp) is the sum of their degrees. We must now introduce in the action large 
N suppressions measured by n(S-^), 

S(T,f) = T-f + y £ j t i N L -^=^ n( - B *' )s > Bi(T,f). (26) 
iei 

It obviously reduces to (TT5|) when n(Bf') = 1 for all I and i S I. The partition function still has the expression (|19l) . 
The amplitude of a connected Feynman graph is 

A(G) = N^t=i •f°e-( d - i )P+ db -E ie r b i Ef=i "(sf ') TT t 6 tj ( 27 ) 

iei 

The formula (|2Tj) still applies to the counting of faces in the subgraphs G">, but it must be used independently for 
each connected component of Indeed, since the subgraphs B^' can have several connected components, G^ 

also typically has several connected components, say n(G^). The latter are (A; + l)-colored graphs and the degree 
w(G^) is the sum of their degrees. By summing the formula (|2"Tj) for each connected component in G®, we get 

J2 /„„ = (A, - 1) P + A, n(G®) (A, - 1) £ 6, n(5«) - ^^ (G* ) + ^ £ 6, u;(5«). (28) 
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When summing over I = 1, . . . , L, the contribution (A/ — i)p becomes (d — L)p and exactly cancels the opposite term 
in (|27[) . The various degrees are positive numbers and their contribution to the scaling is bounded as can be seen by 
applying ([22j) to each connected component in G?W. Therefore we just have to study the terms involving the numbers 
of connected components n(G^),n(B^). The numbers of connected components n{B\ ) are fixed by the choice of 
bubbles, but it is possible to build graphs G with various n(G^). Since this comes with a positive sign in (f28|) . graphs 
with many connected components can be dangerous and might lead to arbitrary high scaling. We therefore need a 
bound on n (G«). 

The fact that G is connected graph can be expressed as follows. If we forget the structure of the bubbles Bi used to 
create G and simply draw them as vertices, then G is a closed connected graph on these vertices, meaning that each 
of them has at least two lines of color which connect it to others. The other lines of color (there are p lines of color 

in total) connect some T to some T which are in the same bubble Bi. Clearly we cannot form n(B^) connected 
components in any Bi that way, because G would then not be connected. We therefore have an optimization problem 
with a constraint (that G is connected). Let us prove the following 

Proposition 1. The maximal number of connected components in the closed, (A; + l)-colored graph G"> is 

n max (G«) = 1 +J2b l {n(B? } ) - l). (29) 



Proof. First we prove that 

»w(G«) > 1 + £>(n(flf ) - 1), (30) 
iei 



by exhibiting a graph G^ with precisely 1 + h( n (B^) — 1) connected components. In each C , we 
choose one connected component and a black and a white vertex in it. Then we form a loop going through all of them. 
This gives one connected component which ensures that G is connected. In each bubble, we are left with exactly 
n(B^) — 1 connected components. They all can become connected components of G^ by putting the lines of color 

between Ts and Ts in the same connected component. After summing over the bubbles, this gives J3iei bi(n(B®) — 1) 
connected components in addition to the one which makes G connected. 

We now use an induction on the number of bubbles in G. The formula of the Proposition [T] trivially holds for one 
bubble (bj = except 6j = 1 for one i El), as n(B^) connected parts can be obtained. Let us then assume that 
the formula holds for any graph built from b bubbles, and consider G to be a graph on b + 1 bubbles, say bj bubbles 
of types Bj, for all j ^ i S / and bi + 1 bubbles of the type Bi for a given i € I. We isolate one of the bubbles Bi, 
which we simply denote B, and distinguish different contributions to n(G^). In the slice I, there are some connected 
components with lines of color only connecting Ts and Ts of B. Denote their number n R . There are some connected 
components which do not go through B, say n^ R of them. Finally, some connected components are shared by B and 
other bubbles, say rS RR . We have n(G^) = + n R + n RR . 

Let us focus on the contributions to n B ^ R . Each of these connected parts goes along some lines of color between 
B and other bubbles. They must actually go along an even number of such lines. Moreover each line of color 

belongs to exactly one connected part. Let us isolate one connected part contributing to n RR . It is a connected 
(A; + l)-colorcd graph. Therefore, cutting the lines of color which connect it to B and reconnecting them in any 
way which disconnects the part included in B from the other bubbles leads to a graph with two connected parts. 
Repeating this process for all connected parts contributing to n RR , all the lines of color which connect B to the 
other bubbles are cut and we see that the number of connected components of G^ is strictly less than the number of 
connected components which can be obtained from B alone plus the number of connected components which can be 
obtained from the b remaining bubbles. Therefore, the formula (|29l) can be used on the b remaining bubbles to get 

n(G«) <n(sf ) ) + l + 6 l (n(Bf ) - l) + ]T b 3 (n{Bf)-\), 

3 ^ eI (31) 
< 2+ (6, + l)(n(sf ) - 1) + E b M B ? ) ~ !)■ 



This states that n(G^) is less or equal to the right hand side of the formula (|2"!)1) with 6j + 1 bubbles of the type Bi 
Combining with (|HTJ)) on bi + 1 bubbles Bi concludes the induction. 
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It comes that 



iei 1=1 



1 jr[A l n(GW)-{A l -l) y £ l b i n(Bj: t >) 
1=1 iei 

which means that the scaling of A{G) with AT is bounded as follows 

\A{G)\ = iTTt- 4 ) ArSf=i[ A '(™(G (0 )-E, e ^Un(sf ) )-l))- T2 ^^(G( ! ') + T3 ^E !e/ b^(B <0 )] 



(32) 



iei 



< ITT^I ^ E ^ i [ T ^ TTT " (G<!>)+Ts ^ TTEie/b, " ( ' Bil!)) ] 



(33) 



The consequence is that the scaling we have chosen in the action restores a balance so that graphs scale uniformly 
(depending on the degrees of subgraphs) , exactly like in the previous section when each G"> has the maximal number 
of connected components and they are slightly suppressed if they have less. 



C. Scaling of observables 

The observables are the expectation values of bubble invariants, 

(B(T, T)) = e F J [dTdT] B(T,f) exp -N d - L S(T,T). (34) 

They admit an expansion over connected Feynman graphs which contain the bubble B as a marked sub-graph. To find 
the scaling of the expectation value, it is sufficient to use the formula we have derived above for an arbitrary graph 
and remove the contribution that would come from the insertion of B as a bubble of the action, i.e. N~ d+ ^'=i n ( Bi \ 
Hence 

(B(T,f)) = N^MB") ^ H£ jj^ 

connected graphs G i£/ 
with d + 1 colors 
containing B 

x N T.tl[^{^)-l-Y. ieI bdn{B^)-l))- J ^u 1 (G^ (35) 

Notice that in the second line the total exponent of AT is negative or at most zero, due to (|2"9"]l . so that at most 
(B) ~ N^'=i ™( s< \ In particular, for the full 2-point function (corresponding to the bubble with exactly two 
vertices), (T ■ T) ~ N L . 



III. THE LARGE N LIMITS 

In the large N limit, the free energy is dominated by graphs G whose subgraphs G"> have vanishing degrees 
u>(G"') = for all I = 1, . . . , L. There are two cases to be distinguished. 

• For a slice of two colors, A; = 2, G"> is a 3-colorcd graph and its degree uj(G^) is the sum of the genera of 
each connected components. Therefore, in the large limit, the connected components of G^ must be planar 
[![. That gives no restriction on the bubbles {Bi}, except that there must exist graphs G whose subgraphs G^ 
for all I are planar at the same time. 

• For slices with at least three colors, i.e. A; > 3, the connected components of G^' must be (A; + l)-colored 
melonic graphs. This can only be true when oj(B^ ) = 0, i.e. the bubbles must be melonic on the colors 
di-i + l,...,di. 

This framework therefore allows more bubbles Bi to contribute in the large AT limit than the scaling of the Section 
U where only melonic bubbles on the colors 1 to d are relevant. For example, at d = 4, with L = 2 and with the color 
decomposition D\ = {1,2}, D% = {3,4}, any bubble can be added to the action and will a priori contribute non- 
trivially, as far as the sub-graphs can be planar at the same time. At d = 5 with L = 2, and D\ = {1, 2}, D2 = 
{3,4,5}, G 1 - 1 ' must be planar and G^ melonic. This implies that all bubbles which arc melonic on the colors 3,4,5 
can a priori contribute (see below why some may not) and in particular the bubbles which are melonic on the five 
colors always contribute. 
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Figure 5. Here d = 6. These two bubbles are melonic on the colors 1,2,3 and on the colors 4,5,6. However, on the left, the 
canonical pairs of vertices for the colors 1,2,3 are (111,1)2) and (113, Va), while there are (1)1,1)4) and (1)2,^3) for the colors 4,5,6. 
Therefore there is no way to add the color while preserving melonicity on the colors 1,2,3 and on the colors 4,5,6 at the same 
time, which makes this bubble irrelevant at large N. On the right of the figure, the canonical pairs coincide for the subgraphs 
with colors 1,2,3 and 4,5,6 (there are (v,v) and (v',v')), so that this bubble is relevant and has a non-vanishing expectation 
value at leading order. 




Figure 6. This is a Gaussian contraction (3 dashed lines of color 0) on a 5-colored bubble (with solid lines). The bubble is 
melonic on the colors 1,2,3 and the contraction is melonic on the colors 0,1,2,3. However, it is not planar on the colors 0,4,5. 
To see it, we count the number of faces: /04 = /05 = /15 = 1, hence the genus is g = 1. The expectation value of this bubble 
has a vanishing contribution at leading order in the Gaussian ensemble with D\ = {1, 2, 3}, D2 = {4, 5}. 



A. Extension of the Universality theorem 



The key point of this universality theorem is that for a given melonic bubble there is only one way to add propagators 
and bubbles from the set {Bi} so as to get melonic graphs. Vertices in a melonic bubble come in pairs (v,v), such 
that the lines connecting v to v carry melonic 2-point insertions. The graphs contributing to the expectation value 
of the bubble have one additional color (the color 0) and melonic contributions are all such that the lines of color 
connect the vertices v and v of the canonical pairs with arbitrary 2-point (melonic) insertions, as shown in the 
Figure |U Therefore, the full contribution is obtained by inserting the full 2-point function of the model between each 
canonical pair of vertices of the bubble. In other words, the expectation value of a bubble factorizes into 2-point 
functions and the models become Gaussian at large N. 

This reasoning holds in our new models on the color slice I when the slice has at least three colors. All graphs G^ 1 ' 
contributing to the expectation value of a bubble B must connect the vertices of the canonical pairs of B^ l \ And it 
must be so for all slices I = 1, . . . , L. This requires that the canonical pairs of vertices are the same for all subgraphs 
in the slices with A; > 3. If this does not hold, the large N limit of the expectation value is zero, and this bubble 
is irrelevant in the action. Examples are given in the Figure O 

If there are in addition some slices with A; = 2, it is also necessary to check that the single melonic contraction 
yields a planar contribution for G^ . It is not always the shown in the Figure [5] 

Therefore, the leading order expectation value (B) ~ ]\[^2i=i n ( B< ) is either zero or Gaussian. Assuming B is a 
connected bubble with 2p vertices with non- vanishing leading order, 



(B(T,T)} =ivSf= 1 «(- B(0 ) [U({U})] P , (36) 
where U = (T ■ T)/N L is the full 2-point function, which depends on the coupling constants {ti}. For instance at 
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which explains why, after re-scaling the bubble on the left by 1/N, these two bubbles can equally contribute to the 
large N free energy, while with the usual scaling, L = 1, the second bubble would not (it is not mclonic on the five 
colors). 

The full 2-point function U({ti}) is obtained explicitly thanks to a Schwingcr-Dyson equation. We start with the 
trivial identity 

E / i dTd ^ ar^— ( r «i-«- cxp-7V d - L S(T,T)) = o. (38) 

ai,...,a d ai---a d 

Then we write explicitly the action of the derivative, which yields 

N d _ N d-L (r.f)-J2ti Pi N d ~^=i (Bi(T, T)) = 0. (39) 

iei 

Notice that all terms have the same scaling N d . Using the factorization (f31)f property, a polynomial equation on U is 
found 

l-U-J2UPiU^ = 0, (40) 

iei 

which is the same equation as ([M]) . This may be surprising because the bubbles Bi are not necessarily melonic on their 
d colors. However, with our scaling, the factorization property (|36j) holds as soon as one subgraph B^ is melonic. 

The result of this section can also be derived by a detailed investigation of the full set of Schwingcr-Dyson equations, 
following the method of |l3| ■ We have checked it explicitly, but do not reproduce it as it is quite long. 



B. Non-Gaussian models at large N 



We have been able to extend the scope of the universality theorem for large random tensors when a color slice A; 
contains three or more colors. But when d is even, it is possible to use a scaling such that all color slices have exactly 
two colors, in which case the theorem does not apply. While we have not solve these models at the end of the day, we 
can prove here that they are not Gaussian at large N . This is the first example of a genuine random tensor model 
which does not satisfy the universality. 

Let d be even, L = d/2 and A; = {21 — 1, 21} for / = 1, . . . , d/2. For a bubble B, n(B^) is the number of loops 
with alternating colors 21 — 1, 21 in B. The free energy F is given by 

cxp -F= f [dTdf] cxp -N d/2 (r-f + Y^U N d/2 -^'=" ™ (s '' >) B t {T, T)) . (41) 

If Bi has a single loop in each color slice, then it receives no extra scaling in the action but the one in front of the 
action. Let us focus on that case to simplify a bit the scalings. The Feynman graphs of that model are (d+ l)-colored 
graphs and the subgraphs are 3-colored graphs, with colors (0,2/ — 1,21). Their degrees are simply the genera 
g( l \G) = (2 — f(G( 1 ') + p)/2, where f(G^) is the number of faces of G^ l K The amplitude of a connected Feynman 
graph in that model is 



(42) 
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(a) This bubble can 
be recast as a 
matrix invariant. 
The non- Gaussian 
behavior is thus 
well-known. 




(b)This is a genuine 

tensor invariant 
because the loop of 
colors 34 is twisted 
with respect to the 
loop with colors 12. 



Figure 7. Two 4-colored bubbles which lead to non-Gaussian behaviors. 




(a) / 03 = 3, (b) / 03 = 2, (c) / 03 = 2, (d) / 03 = 2, (e) / 03 = 1, 

/04 = 1, 3034 = /04 = 2, 9034 = fo4 = 2, 9034 = /o4 = 2, 9034 = fo4 = 1, 9034 = 1 



Figure 8. The graphs represent Gaussian contractions of a bubble with propagators (dashed lines of color 0). The sub-graphs 
with colors 012 are always planar. Below each graph we write the number of faces of colors 03 and 04. The genus of the 
sub-graphs with colors 034 are found by go34 = (2 — /03 — /04 — 1 + 3)/2. It shows that the contraction on the far right is not 
planar on the colors 034. 



In other words, one has the scaling of a matrix model for each subgraphs . This implies that at large N, there 
is no melonicity involved anymore, but only planar graphs. Since it is well-known that planar graphs typically lead 
to non-Gaussian large N behaviors, this class of models looks promising to escape the universality of large random 
tensors. 

We now show that this is true with an example. We set d = 4. A simple example of bubble is the one in the Figure 



7(a) Notice however that it can be recast as a matrix invariant of order 6 for a matrix of size N 2 x TV 2 . The matrix 
is Mab = T aia2a30ii with A = (01,03) and B = (02,04) and the invariant is then tr(MA-/^) 3 . Then our tensor model 
is just a matrix model, which can be solved using matrix model techniques. Consider instead the bubble with one 
loop of color (12) and one loop of color (34) shown in the Figure [7(b)| This bubble is really a tensor invariant and 
cannot be recast as a matrix invariant. 

A planar contraction on this bubble is shown in the Figure |8(a)| which scales like N 4 as expected. Another 
planar contribution is shown in the Figure [9] It has 2 bubbles, 6 lines of color and 12 faces, hence scales like 
jy2x2-6x2+i2 _ ^y4_ Moreover, it clearly cannot be displayed in terms of two-point insertions on one of the two 
bubbles but it is instead a genuine 6-point contribution to the expectation value of the bubble. Therefore this 
contribution is not Gaussian. 

Finally, we insist on the following point. The dominant graphs in the models presented in this section must have 
their subgraphs with colors 0,1,2 planar, with colors 0,3,4 planar and so on. In the model built with the bubble from 
the Figure 7(a)| the subgraphs with colors 0,3,4 are planar as soon as as the subgraphs with colors 0,1,2 are planar 
(because they simply are the same). In particular, the expectation value of this bubble in the Gaussian ensemble 
of covariance A is the same (after rescaling by A -2 ) as the expectation value of tr(MM'l') 3 (rescaled by A -1 ), i.e. 
A 3 x 5, where 5 is the Catalan number C3 counting the number of planar contractions However, with the bubble 
of the Figure [7(b)| the subgraphs with colors 0,1,2 and with colors 0,3,4 are different and may not be planar at the 
same time. For instance, in the Gaussian ensemble of covariance A, the expectation value of this bubble (rescaled by 



A 2 ) is A 3 x 4. The Figures 8(a)[ |8(b)[ |8(c) [ |8(d)[ |8(c) indeed shows that only four of the five planar contractions for 



the subgraph with colors 0,1,2 are also planar for the subgraph with colors 0,3,4. 



13 




Figure 9. This is a two bubble graph which is planar on the colors 012 and 034, and therefore scales like the Gaussian planar 
contractions in the Figure [8] But it is a non-Gaussian contribution to the free energy. 



CONCLUSION 



In this paper, we have generalized Gurau's 1/N expansion to a family of well-defined expansions. There are more 
available expansions as the dimension d increases (the new scalings we are proposing appear at d = 4). The key 
idea is to slice the set of indices of the random tensor into groups of two or more indices and associate to a colored 
graph the degrees of its sub- graphs determined by the slicing, as described in the Section Ql] This re-organizes Gurau's 
expansion, enhancing the contribution of previously suppressed graphs, in such a way that for instance more than 
the melonic sector contributes at large N. However, as soon as there is a slice with at least three colors, the large N 
behavior can only be Gaussian, because the corresponding sub-graphs must be melonic. This result is an unexpected 
extension of the universality theorem of [l2|. 

One motivation of this work was to identify scalings leading to non-Gaussian large N limits, thereby stepping out 
of the universality. This is indeed possible for even d, by looking at the genus of ribbon sub-graphs carrying only 
some subsets of colors. Therefore, it seems that such a scaling is closer to matrix models than to the large N limit of 
random tensors known so far. The next step is thus to find the large N solution (and eventually the solution at all 
orders), describe its continuum limits and critical exponents. While we expect them to be different of the universality 
classes already found in tensor models, it is not clear whether it will coincide with the critical behaviors of matrix 
models or be completely new. 

A question that has not been discussed in the main text of the paper is the topology of the graphs. Indeed, the 
quantum gravity motivation in random tensors is that the bubbles are dual to (d — l)-dimensional pseudo-manifolds 
and the Feynman graphs in the expansion of the free energy are dual to d-dimensional pseudo-manifolds. It is 
known that the melonic graphs are dual to triangulated spheres, generalizing the relation between planar graphs and 
discretizations of the sphere in two dimensions. However it is also known that not all colored triangulations of the 
sphere are melonic. Quite evidently, our new large N limits pick up contributions from spheres which are not melonic, 
which is necessary to get more realistic quantum gravity models. But in addition, our large N limits may contain 
non-spherical topologies (e.g. when there is a slice with two colors only). 

Finally, we would like to mention a more speculative idea. It has been so far observed that although the degree 
reduces to the genus at d = 2, the combinatorics of the dominant graphs is quite different between matrix and tensors 
for d > 3. The physical consequences are the universality theorem for d > 3, different statistical behaviors and 
hence different critical exponents in the continuum limit. However, the present paper shows that tensor models admit 
matrix-like scalings (i.e. with the genus of ribbon sub-graphs). Therefore, it is natural to ask whether interpolating 
between matrix-like and tensor-like behaviors is possible. While appealing, this idea seems difficult to realize. It 
would certainly necessitate some analytical continuations on d, L (the number of color slices) , and on the number of 
colors in the slices. The difficulty is that the difference between d — 2 and d > 3 only appears in the formula for the 
degree © and it involves faces, i.e. closed subgraphs labeled by pairs of colors. It is unclear how to go to non-integer 
numbers of colors from that view. 

If this question can be addressed, we think it could be through the Schwingcr-Dyson equations. Indeed, the set of 
Schwinger-Dyson equations found in jl8| is formally the same as for the new scalings. They involve the same set of 
operations on bubbles. Only the powers of N are modified so as to change the balance between the different orders 
in the 1/N expansion of the bubble expectation values. We plan to study this set of equations for the non-Gaussian 
scaling found in this paper and compare with the large N equations obtained and solved in [ill ]. Furthermore, all 
solutions to tensor models are certainly based on this set of equations with different balances, and this is where an 
interpolation between matrix-like and tensor-like behaviors could be found. 
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